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INTRODUCTION 


Uncertainties  about  turbulent  combustion  in  hydrogen-air  systems  have  an  impact 
on  our  abilities  to  achieve  successful  designs  of  supersonic  combustion  devices  for  the 
National  Aerospace  Plane,  for  example.  Prior  to  initiation  of  the  present  research,  it  was 
not  known  whether  the  combustion  was  more  likely  to  occur  in  reaction-sheet  or 
distributed-reaction  regimes,  and  common  practice  was  to  assume  the  latter  and  employ 
laminar  chemical  kinetics  in  routines  based  on  computational  fluid  mechanics.  As  will  be 
seen  below,  estimates  suggest  that  this  is  a  poor  procedure  since  the  reaction-sheet  regime 
is  more  prevalent  than  the  distributed-reaction  regime.  The  present  research  was  initiated 
because  it  was  felt  desirable  to  try  to  determine  the  regimes  more  firmly  and  to  develop 
improved  methods  for  accounting  for  the  chemical  kinetics  of  the  combustion. 

RESEARCH  OBJECTIVES 

The  objective  of  this  research  is  to  improve  understanding  of  the  chemical 
kinetics  and  fluid  dynamics  of  turbulent  combustion  in  high-speed  flows.  Supersonic 
combustion  in  hydrogen-air  mixtures  is  being  addressed  by  theoretical  approaches  that 
distinguish  between  reaction-sheet  and  distributed-reaction  regimes.  The  work  seeks  to 
identify  effects  of  compressibility  in  turbulent  combustion,  methods  for  including 
compressibility  in  theoretical  analyses,  and  reduced  chemical-kinetic  mechanisms 
appropriate  for  supersonic  combustion.  The  results  may  help  to  enhance  capabilities  of 
reasonable  computations  of  high-speed  turbulent  reacting  flows. 

ACCOMPLISHMENTS 

Regimes  of  Turbulent  Combustion 

The  work  on  the  present  project  has  been  addressing  flamelet  regimes  rather  than 
distributed-reaction  regimes  of  turbulent  combustion.  The  motivation  for  focusing  on 
flamelets  stems  from  the  initial  evaluation,  under  the  previous  grant,  indicating  that  in  the 
application  to  aerospace  propulsion,  the  hydrogen-air  system  lies  in  the  reaction-sheet 
regime.  Those  initial  estimates  were  made  prior  to  development  of  the  extensive  data 
base  on  hydrogen-air  diffusion  flames,  from  the  previous  grant.  Because  of  this,  it 
seemed  reasonable  at  the  outset  to  reevaluate  the  combustion  regimes,  on  the  basis  of  the 
new  information  that  has  now  been  accumulated.  The  results  of  this  reevaluation  reaffirm 
that  reaction-sheet  flamelet  regimes  are  most  relevant.  In  comparison  with  the  original 
estimates,  there  has  been  a  narrowing  of  the  most  probable  range  and  a  displacement 
towards  the  distributed-reaction  regime  as  a  consequence  of  improved  chemical-kinetic 
and  combustion  information.  Therefore,  attention  continues  to  be  placed  on  flamelets. 

Ignition  in  Hvdrogen-Air  Mixtures 

Previous  work  addressed  steady  burning  and  extinction  in  hydrogen-air  systems 
but  did  not  consider  the  ignition  branch  of  the  diffusion-flame  S  curve.  Autoignition  of 
hydrogen  with  hot  air  is  of  fundamental  interest  in  high-speed  propulsion;  for  instance,  in 
a  supersonic  combustor  the  fuel  jet  is  expected  to  be  injected  into  and  ignited  by  a  stream 
of  hot  air  heated  by  shocks.  Typical  studies  of  ignition  focus  on  the  chemistry  with  little 
or  no  attention  being  paid  to  the  fluid  dynamics.  Complex  interactions  which  are  likely 


to  play  art  important  role  in  the  ignition  processes  in  real  flow  are  thus  ignored.  To  help 
to  understand  these  interactions,  current  attention  is  focused  on  the  ignition  of  a  cold  jet 
of  hydrogen  impinging  against  an  opposed  jet  of  heated  air. 

The  problem  is  reduced  to  a  two-point  boundary-value  problem  for  ordinary 
differential  equations.  The  numerical  solution  to  the  discretized  equations,  normally 
obtained  employing  New'ton's  method  or  a  variant  thereof,  encounters  difficulties  in  the 
neighborhood  of  tumirg  points  where  the  derivative  of  the  solution  with  respect  to  a 
parameter  of  the  problem  becomes  infinite  (for  example,  at  the  extinction  condition).  To 
circumvent  these  difficulties,  a  modification  of  the  basic  algorithm  is  introduced  in  the 
present  work,  employing  a  path-tracing  continuation  method  of  Smooke,  in  which  points 
on  the  arc  of  solutions  are  computed.  This  method  has  now  succeeded  in  generating  the 
ignition  branch. 

Computations  began  with  the  21 -step  mechanism  which  was  used  for  extinction 
studies  reported  earlier.  Numerical  integrations  were  also  performed  with  a  different 
chemical  scheme  proposed  by  Smooke  and  his  coworkers.  For  dilute  hydrogen  flames  at 
one  atmosphere  pressure,  a  large  discrepancy  in  the  extinction  and  ignition  strain  rates 
was  observed  with  increasing  feed  temperature  of  the  air  stream.  Furthermore,  thermal 
diffusion  was  seen  to  play  a  surprisingly  important  role  in  diluted  flames.  In  view  of 
these  unexpected  preliminary  results,  a  critical  assessment  of  the  reaction  rates  was 
deemed  warranted.  Af.er  an  extensive  survey  of  the  literature  and  discussions  with  other 
researchers,  especially  Yetter,  the  21-step  chemical-kinetic  scheme  with  associated  rate 
parameters  shown  in  Table  1  was  selected  as  representing  the  best  information  currently 
available. 

To  test  the  updated  chemistry,  flame  speeds  for  hydrogen-air,  wet  CO  and 
methane-air  mixtures  were  calculated.  The  agreement  between  theory  and  experiment  is 
excellent  over  a  wide  range  of  initial  conditions  for  the  hydrogen-air  flame.  A  flame 
speed  of  38.6  cm/s  for  stoichiometric  methane-air  flames  was  obtained,  in  good 
agreement  with  the  experimental  value  of  38-40  cm/s.  Of  particular  interest  in  the  new' 
chemical  scheme  is  the  rate  of  the  chain-branching  reaction  H  +  O2  -»  OH  +  O,  which  is 
based  on  the  recent  Stanford  shock-tube  experiments.  The  new  rate  provides  a  slightly 
larger  radical  pool,  and  unlike  the  rate  recommended  by  Baulch  et  al.  in  their  recent  1992 
review,  exhibits  the  correct  direction  of  the  temperature  dependence  for  the  reverse 
reaction  OH  +  O  H  +  O2  in  the  temperature  range  of  500-3000  K. 

Calculations  to  characterize  the  ignition  of  hydrogen-air  mixtures  over  a  wide 
range  of  conditions  of  interest  in  high-speed  propulsion  devices  are  now  being  performed 
with  the  updated  chemical  scheme. 

Structures  of  Counterflow  Diffusion  Flames  with  Small  Stoichiometric  Mixture  Fractions 

Asymptotic  methods  were  used  to  evaluate  a  characteristic  diffusion  time  for 
counterflow  diffusion  flames,  a  quantity  of  utmost  importance  in  analyses  of  flame 
structure  and  extinction.  Attention  was  restricted  to  the  important  limit  of  small  values  of 
the  stoichiometric  mixture  fraction,  which  resulted  in  large  values  of  the  stream  function 
at  the  reaction  sheet,  a  quantity  which  was  then  employed  as  a  large  expansion  parameter. 
Three  distinct  layers  were  identified  inside  the  viscous  layer:  (1)  a  convective-diffusive 


layer  on  the  oxidizer  side  of  the  reaction  sheet,  (2)  a  rotational  inviscid  layer  on  the  fuel 
side  of  the  reaction  sheet,  and  (3)  a  fuel-product  mixing  layer  in  the  vicinity  of  the 
stagnation  plane;  all  three  of  these  layers  were  analyzed  and  matched.  The  matching 
produced  analytic  expressions  for  the  stream  function  and  mixture  fraction  at  the  reaction 
sheet  and  corrected  formulas  for  the  characteristic  diffusion  time  previously  derived  on 
the  basis  of  constant-density  or  nonreactive-mixing  assumptions.  Comparison  between 
analytical  and  numerical  results  demonstrate  that  the  analytical  results  obtained  from  the 
present  analysis  provide  accurate  expressions  for  the  characteristic  diffusion  times  and 
stoichiometric  scalar  dissipation  rates.  The  results  provide  improved  predictions  of  the 
diffusion-flame  structure  and  extinction. 

Nonpremixed  Flames  in  Stagnating  Turbulence 

It  is  becoming  widely  recognized  that  the  study  of  flames  in  stagnating  turbulence 
nrovidcs  a  powerful  vehicle  for  examining  fundamental  issues  in  turbulent  ^nii"Kt:nn 
For  the  present  discussion  the  flow  of  interest  consists  of  two  coaxial,  axisymmetric  and 
closely  spaced  jets  within  which  a  turbulence-generating  grid  is  installed.  By  changing 
the  mean  velocity  from  each  jet,  the  entire  range  of  chemical  behavior  from  equilibrium 
to  extinction  can  be  covered  in  one  experimental  setup.  The  initial  experimental 
investigations  of  these  flames  were  confined  to  premixed  systems,  i.e.  to  a  turbulent 
reactant  stream  either  impinging  on  a  wail  or  opposing  another  reactant  stream,  but 
recently  preliminary  results  of  an  experiment  involving  a  flame  in  a  turbulent  methane 
stream  opposing  a  turbulent  airstream  have  been  reported.  Since  the  only  theory  of 
flames  in  stagnating  turbulence  is  concerned  with  premixed  systems,  an  analysis  of  the 
nonpremixed  case  has  been  started  under  the  present  program. 

The  initial  results  are  based  on  the  thin-flame  approximation  and  on  chemical 
equilibrium,  with  matching  to  outer  zones  imposed  on  each  side  of  the  flame.  This 
approximation  involves  two  nonreactive  streams  impinging  on  the  opposite  sides  of  a 
turbulent  flame  whose  thickness  is  small  compered  with  the  jet  spacing.  The  resultant 
calculated  temperature  distribution  exhibits  a  peak  value  and  general  shape  in  excellent 
agreement  with  experiments,  but  the  predicted  thickness  in  physical  space  is  an  order  of 
magnitude  too  small.  In  brief,  the  thin-flame  approximation  is  not  applicable  to  the 
experimental  conditions.  The  combination  of  heat  release  and  of  turbulent  length  scale 
preclude  quantitative  comparison  of  theory  with  these  experiments.  Accordingly,  an 
alternative  analysis  has  been  developed,  describing  in  one  set  of  equation  the  entire  flow 
from  one  jet  exit  plane  to  the  other.  The  numerical  treatment  of  this  new  set  of  equation 
is  underway.  Initial  application  will  be  to  the  methane-air  system,  to  permit  comparison 
with  available  experimental  data;  if  the  agreement  is  satisfactory,  with  or  without 
refinements  to  the  theory,  the  analysis  will  then  be  applied  to  the  hydrogen-air  system. 
Additional  studies  along  these  lines  involve  the  effects  of  finite  chemical  rates,  so  that 
extinction  and  ignition  can  be  analyzed.  These  studies  will  clearly  relate  to  the  flamelet 
analyses  described  earlier. 
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Table  1:  Specific  Reaction-Rate  Constants  for  the  Hydrogen  Oxygen  Reaction  System  Adopted  in  the  Present 
Study. 


No. 

Reactions 

A1 

n1 

E1 

1 

H  +  02  s=s  OH  +  O 

3.52x  1016 

-0.7 

17070 

2 

H2  +  O  ^  OH  +  H 

5.06xl04 

2.67 

6290 

3 

OH  4*  OH  ^  H2O  +  O 

Jfc  =  5.46  x  10u  exp  (0.00149r) 

4 

H2  +  OH  H20  +  H 

1.17xl09 

1.3 

3626 

52 

H  +  02  +  M  —  H02  +  M 

6.76xl019 

-1.42 

0 

6 

H  +  H02  —  OH  +  OH 

1.70xl014 

0.0 

874 

7 

H  4*  HO2  — *  H2  4*  O2 

4.28xl013 

0.0 

1411 

8 

OH  +  H02  — *  H20  +  02 

2.S9xl013 

0.0 

497 

93 

H  +  H  +  M  — ♦  H2  +  M 

l.SOxlO18 

-1.0 

0 

102 

H  +  OH  +  M  -*  H20  +  M 

2.20x  1022 

-2.0 

0 

11 

HO2  4"  HO2  — *  H2O2  4-  O2 

3.02X1012 

0.0 

1390 

124 

H2O2  +  M  —  OH  +  OH  +  M 

1.20xl017 

0.0 

45500 

13 

H202  +  OH  H20  +  HO2 

7.08xl012 

0.0 

1430 

14 

O  +  H02  ^  OH  +  02 

2.00xl013 

0.0 

0 

15 

H  +  HOz  =*  O  +  H20 

3-lOxlO13 

0.0 

1720 

162 

H  +  O  +  M  —  OH  4-  M 

6.20x  1016 

-0.6 

0 

172 

O  4  O  4  M  ^  O2  4*  M 

6.17xl015 

-0.5 

0 

18 

H2O2  +  H  ~  H2O  +  OH 

l.OOxlO13 

0.0 

3590 

19 

H202  +  H  HO2  +  H2 

4.79xl013 

0.0 

795C 

20 

O  4  OH  +  M  —  H02  +  M 

l.OOxlO16 

0.0 

0 

21 

H2  +  02  —  OH  +  OH 

1.70xl013 

0.0 

47780 

“  Units:  mol/cm3,  s-1,  K,  cai/mol;  rates  for  reverse  steps  obtained  from  JANAF  thermo¬ 
chemical  equilibrium  data. 

1  Chaperon  Efficiencies:  H2:  2.5,  H20:  12.0,  O2:  1.0  and  N2:  1.0 
e  Chaperon  Efficiencies:  H2:  1,  H2O:  6.5,  O2.'  0.4  and  N2:  0.4 
d  Chaperon  Efficiencies:  H2:  2.5,  H2O:  15.0,  02:  1.0  and  N2:  1.0 
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